Universality of the Crossing Probability for the Potts Model for q = 1,2,3,4 
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The universality of the crossing probability tths of a system to percolate only in the horizontal 

direction was investigated numerically by a cluster Monte-Carlo algorithm for the g-state Potts 

model for q = 2, 3, 4 and for percolation q — 1. We check the percolation through Fortuin-Kasteleyn 

clusters near the critical point on the square lattice by using representation of the Potts model as 

the correlated site-bond percolation model. It was shown that probability of a system to percolate 

only in the horizontal direction tths has universal form tths = A(q)Q (z) for q = 1, 2, 3, 4 as a function 

r i i<(«) 

of the scaling variable z = \b(q)L" i .i'> (p — pc(q, L))\ . Here, p = 1 — exp(— /3) is the probability 

of a bond to be closed, A(q) is the nonuniversal crossing amplitude, b(q) is the nonuniversal metric 
factor, i/(q) is the correlation length index, ((q) is the additional scaling index. The universal 
function Q(x) ~ exp(— \z\). Nonuniversal scaling factors were found numerically. 

PACS numbers: 64.60.Ak, 05.10.Ln, 05.70.Jk 



I. INTRODUCTION 

The concept of the universality and scaling relations [l| 
is general concept of the modern phase transition theory. 
The main point of the scaling theory is that in the vicinity 
of the critical point for a system of linear size L, the crit- 
ical behavior of thermodynamical quantities can be ex- 
pressed as a universal function of two variables: reduced 
temperature r = T ~ Tc and external field h. The finite- 
size scaling of thermodynamical functions of spin models 
was studied theoretically and numerically |2|, bJ U, |5j . In 
Ref. 0, Privman and Fisher proposed the idea of the 
universal finite-size scaling with nonuniversal metric fac- 
tors. For example, for the free-energy density of system 
size L and dimension d 

f{T, h; L) = L- d Y (c 1 TLi,C 2 hL ( - l3+ ' l) /^ (1) 

where v is the scaling index for correlation length, (3 
is the scaling index of magnetization and 7 is the scal- 
ing index of magnetic susceptibility. For systems of dif- 
ferent boundary conditions, aspect ratios and geome- 
tries (square, honeycomb, triangle), the scaling function 
Y(x,y) is universal and only metric factors C\, C2 are 
system dependent. Scaling properties of thermodynam- 
ical functions of the Potts model were investigated in 
Refs. H@. 

Langlands, Pichet, Pouliot and Saint-Aubin 9] show 
that for site and bond percolation on square, honeycomb 
and triangle lattices with the aspect ratios a, and 
av3/2 respectively, the crossing probability -kk of a sys- 
tem to percolate in the horizontal direction is the univer- 
sal function of a. Hu, Lin show that by choosing a very 
small number of nonuniversal metric factors, all scaled 
data for percolation functions and the number of perco- 
lating clusters on square, honeycomb and triangle lattices 



may fall on the same universal scaling functions 0, ^| . 
Their scaling argument was x — (p ~ p c )L^ where p c is 
the critical point, L is the lattice size, v is the correla- 
tion length index. The scaling of crossing probabilities 
for the three-dimensional percolation was investigated in 
Ref. [12. 

The continuum limit of the crossing probability 717, (p c ) 
was investi gate d by J. Cardy by conformal field meth- 
ods Theanalo gous formula for the crossing 
probability -Kh v = tth — W hs w as found by Watts [T^ |. 
The works of Smirnov [TR Il8| analytically proved that 
the crossing formula holds for the continuum limit of site 
percolation on the triangle lattice 0, 0] . 

The g-state Potts model can be represented as the 
correlated site-bond percolation in terms of Fortuin- 
Kasteleyn (FK) clusters [3 ■ At the critical point of the 
second order phase transition Potts model, FK-clusters 
exhibit the percolation transition. So there is an intrin- 
sic relationship between critical properties of the Potts 
model and percolation properties of FK clusters. The 
universality of the crossing probability for the Ising model 
on rectangle lattices of square, honeycomb and triangle 
geometries was investigated by Langlands et all. |2(J, |21j . 
Hu, Chen and Lin show the universality of the crossing 
probability and a number of percolation clusters for 
the correlated site-bond percolation q = 2 . 

In this paper the the crossing probabilities of FK clus- 
ters is studied numerically. We investigate the universal- 
ity of the crossing probability with respect to a number 
of spin states q of the Potts model. We show numerically 
that the probability of a system to percolate only in hori- 
zontal direction 7r/ ls is an universal function of the scaling 



variable z = 



C(«) 



for q = 1, 2, 3, 4 where 
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b(q)L»(p-p c ) 

p = 1 — exp(— j3) is the probability of bond to be closed, 
/3 = 1/TkB is an inverse temperature, p c = 1 — exp(— j3 c ) 
is the location of critical point in the p scale, b{q) is a 
nonuniversal metric factor and Q(q) is a scaling index, 
described dependence of the form of the crossing proba- 



2 



bility on q. 

We show that for each value q — 1, 2, 3, 4 on the square, 
lattice the index £ is practically does not depend on the 
lattice size. Therefore, by introducing this index £(<;), we 
can lay on the same curve the points for different q in the 
critical region. 



II. CROSSING PROBABILITIES FOR THE 
POTTS MODEL 

In the Potts model, each spin <jj can take one of the 
q different values I, . . . , q and the Hamiltonian is written 
as Ti — — J ^2 5(<Tj, (Tj), where J is the ferromagnetic 

coupling constant, which we set it equal 1. The partition 
function of the Potts model j23J is 

Z = 5>xp(-/3W(u;)) = J2 II K 1 ~P) +P 5 (°i^i)\ 

a a (ij) 

(2) 

where (3 — 1/Tks is the inverse temperature, p — 
1 — exp(— (3) is the probability of bond to be closed 
and 1 — p = exp(— (3) is the probability of bond to be 
open; summation is performed over all spin configura- 
tions a. Sometimes the term in square brackets is ex- 
pressed via v = exp(/3) — 1: [(1 — p) + pS(a{, aj)] = 
exp(— f3) [1 + v5(<Ti, (Tj)], but we write Eq J2J via p to 
emphasize the fact that in correlated site-bond percola- 
tion ^3 , the probability of bond to be closed is p and to 
be open is 1 — p. 

For the square lattice of linear size L with periodic 
boundary conditions the total number of bonds is N — 
2L 2 . Let us define by L the graph of all the edges on 
the lattice. The product over all bonds consists of 
N terms, so the product may be expanded as sum of 
2 N terms. Let us associate to each of these 2 N terms 
subgraph G of graph £, by following rule. Each of 2 N 
terms can be considered as the product of N factors. 
Each factor [(1 — p) + p8(<Ti, corresponds to some 
edge of the graph C. To construct G, we perform 
the following procedure. If this factor for edge is 
equal to 1 — p, we delete edge from subgraph. If 

this factor for edge is equal to p, we leave this 

edge in the subgraph. So we associate to each term 
in sum (|2J subgraph G. Each subgraph G consists of 
6(G) edges (closed bonds) and G(G) connected compo- 
nents. The term in © corresponding to G contains factor 
(1— p) N ~ b ( G )p bt y G \ and delta-functions guarantee equiva- 
lence of spins in each connected component. As a result of 
the summation over all possible spin configurations, the 
contribution of this subgraph G into the partition func- 
tion is q c ( G \l— p) N ~ b ( G )p b ( G \ Therefore, we can replace 
the summation over all spin configurations in Eq. J2J by 
summation over all possible subgraph G on C 

z=Y,<f [G) o--p) N ~ h{G) p b{G) ( 3 ) 

Gee 



We shall keep in mind that this definition is valid even for 
the non-integer value of q. The partition function of the 
Potts model for q — > 1 corresponds to the percolation, 
where p is a probability of bond to be occupied. 

Now we introduce the crossing probability 7r; ls , the 
probability of a system to percolate only in the horizontal 
direction while the percolation in the vertical direction is 
absent. We must distinguish it from the probability of a 
system to percolate in the horizontal direction irrespec- 
tive to percolation in the vertical direction tth . We define 
the indicator function I; ls (G) for each subgraph G in ac- 
cordance with the rule 

I (G) f 1 if C percolates only in horizontal direction 
hs{ J - | o in all others cases 

(4) 

We mean that G percolates only in the horizontal di- 
rection if it contains at least one connected component 
touching left and right sides of the lattice, and there are 
no components joining top and bottom sides. Therefore, 
the crossing probability 717^ (/3; L,q) of a system of size 
L with q possible spin states to percolate only in the 
horizontal direction at [3 can be written as 

*H.(PiL,q) = ± £ hs{G)q C ^{l-p) N ^ G V {G) (5) 
Gee 

Here, p = 1 — exp(— (3). We can introduce the crossing 
probability for vertical direction ir vs by the same way. 
But on the square lattice -Khs = ~k vs and later on, inves- 
tigate only ir hs . 

Let us notice that 7r/ ls (/3; L, q) has maximum near the 
critical point (3 C (L, q) because in the ordering phase f3 > 
p c (L,q), the most probable configurations contain large 
clusters touching top and bottom sides of the lattice, and 
do not contribute into 7r^ s . 



III. NUMERICAL RESULTS 

We use the Wolff cluster algorithm [24| to generate the 
different spin configurations. For each spin configuration 
we generate a bond configuration in accordance with . 
Then, we break the lattice into independent clusters of 
connected sites by using the Hoshen-Kopelman [25| al- 
gorithm. Then, we analyze crossing properties Ih s (G) of 
this configuration. Between checking the crossing we skip 
5 Monte-Carlo steps. For each value of /?, the averaging 
is performed over 10 series each of 10 5 configurations. 
The total number of configurations is 10 6 . Sets of con- 
figurations are used for the estimation of the numerical 
inaccuracy. Quantity Its is an indicator function. It 
means that it takes values or 1 for each configuration. 
Therefore, the resolution of our computations for tt^s is 
10~ 6 . We compute data for q = 1 (percolation), 2, 3, 4 
and lattice sizes L = 16,32,48,64,80,96,112,128. For 
each pair (q,L), we perform computation for 200 values 
of (3 (or p for percolation) in the critical region. 
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There is a very important question, how the scaling 
variable should be chosen? 

For the percolation q = 1 the scaling variable is defined 
as the deviation of the bond concentration from critical 
point r = p — p c . For the Potts model, we can take 



3. r = (1 - exp(-/3)) - (1 - exp(-/3 c ) = p - p c 

The first variant t = T ^ Tc is usual for the investigation 
of thermodynamic quantities 0, |23, the second one is 
widely used for the approximation of magnetic suscepti- 
bility x °f the Ising model in the critical region [2(| , and 
the third variant r = p — p c may be argued by the fact 
that for the site-bond correlated percolation, the proba- 
bility of bond to be occupied is p = 1 — exp(— (3) - see 
equations 0,0- 

Let us evaluate the symmetry of the crossing probabil- 
ity as a function of all these variables by comparing third 
moments of the function 7r/ ls . We assume that iThs must 
be the most symmetric function as a function of " right" 
variable. We calculate four moments of -Khs{x) as a func- 
tions of x = £, t, r by numerical integration in accordance 
with formula uq = J ^hs(x)dx, Ml = 7^ / XTTf ls (x)dx ) 
fik = J(x — fii) k TTi ls (x)dx. Results of computations 
//3 for the lattice size L = 128 are given in Table |U We 

TABLE I: The third moment /i3 of the crossing probability 
■Kha, computed for different variables t, r and r 



q 


2 


3 


4 


t 


-7.1(3) x 1CT 8 


-9.1(2) x 1(T 9 


-8(1) x ict 10 


T 


-2.1(6) x 10" 8 


-2.3(2) x 1(T 9 


-2(77) x ICT 12 


r 


2.5(2) x 10~ 9 


2.8(1) x KT 11 


1.3(5) x ICT 11 



see that for all the cases ^3 is smaller for scaling variable 
r, with one exception q — 4, r, when numerical error is 
approximately thirty times greater than value /Z3 (r) and 
six times greater than Hz(r). For all other lattice sizes 
L = 16,.. .,112 the third moment /X3 calculated by us- 
ing variable r, is smaller than for t or for t. Therefore, 
we will work with the scaling variable p = 1 — exp(— 0), 
the probability of bond to be closed for q = 2, 3, 4. The 
critical point of the q-state Potts model on the infinite 
lattice is [23j p c (q) — Ty^~pr- The bond percolation critical 

point p c {q = 1) = 5 can be obtained from this formula. 
Variable p = 1 — exp(— (3) naturally provides crossover 
from the percolation to the Potts model. We can plot 
the crossing probability TTh s as a function of p for the 
percolation and the Potts model on the same graph. 

In FigE] data for nf ls (p; L, q) are plotted for q = 
1,2,3,4 and L = 32,64,128. We can see shift of the 
critical point p c (q,L) for finite lattice sizes L. We also 



'L=32 




FIG. 1: The crossing probability Hhs(p\ L,q) for q = 1,2,3,4 
and L = 32, 64, 128 

see the change of the width of the function due to size 
scaling. Let us try to identify the shape of 7r^ s by an- 
alyzing the ratio We know that at least for the 
percolation, function 7T/,. S looks like the Gaussian func- 
tion ~ Aexp(— x 2 ), and for the Gaussian function 
we expect (Gaussian) = 3. But, in reality, cross- 
ing probabilities are not Gaussian. The moments ratio 
for function Tr' h (p) (the derivative of crossing probabil- 
ity Tthip) with respect to p), for the percolation model 
was found by Ziff U$(q = 1) ~ 3.20(5) 28] and more 

precisely p{q = 1) ~ 3.176(4) 29] and by Hovy and 
Aharony (q = 1) ~ 3.174(25) 30]. The explanation 
of the non-Gaussian crossing probability shape was first 
given by Berlyard and Wehr [3l| and was also discussed 
in Ref. [33]. Newman and Ziff verify that the tails of 
the crossing probability falls of as exp(— c(p — p c ) 4 ^ 3 ), 
therefore the tail of the distribution is characterized 
by the correlation length exponent [3^, 0]. However, 
as the first approximation we introduce an additional 
scaling index £ and check moment ratios for function 
g(x;C) = Aexp(— x^). For this function g(x;£) the mo- 
ment ratios are 

We calculate moments for several values of £ and put 
these data into the first (£) and second (^(0) rows of 
the Table fill The choice of values £ will be argued later. 
Results of computation of the moments ratio (a) for 

the crossing probability TTh s (p; L,q) for L = 32,128 is 
placed in the fourth and fifth row of the Table [n] We 
see that the moments ratio practically does not depend 
upon the lattice size L. We check this fact for others lat- 
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tice sizes. Moment ratios for g(x\ Q slightly differ from 

TABLE II: The ratio of moments for g(x,0 and 
n h s(p;L,q) 



c 


2 


3/2 


4/3 


6/5 


£f (Q analytically 


3 


3.76195 


4.22219 


4.7434 


q 


1 


2 


3 


4 


^ (q) numerically L = 32 


3.145(6) 


3.871(7) 


4.577(14) 


5.28(2) 


U%(q) numerically L = 128 


3.18(2) 


3.91(8) 


4.56(2) 


5.30(4) 



moment ratios for n^s (this fact will be explained later), 
but we can try to approximate the crossing probability 
by the function g(x, £) and then compare results of ap- 
proximation with numerical data. 

Below, we describe the fitting procedure. As we can 
see from Fig^ there are many nonuniversal scaling fac- 
tors for 717^: the amplitude of the crossing probability 
A(L,q), the finite size critical point p c (L,q), which dif- 
fers from p c (L ~ 00, q) = , the nonuniversal scaling 

factor B(L,q), which provide the finite size scaling of the 
function, and the additional scaling index ((L,q). 

We perform the four-parametric fit of 7r^ s (p; L,q) by 
the function Q) 

F(p; L, q) = A(q, L) exp(- [B(L, q)(p - p c (L, q))f L ' q) ) 

(7) 

We use points log(7r/ ls ) > —9 for this fit and the log 
scale for the ordinate axis. As a result of the approxima- 
tion we get a set of scaling amplitudes A(q,L), nonuni- 
versal metric B{q 1 L) 1 critical points p c (q,L) and scaling 
indices ((q,L). Then, we use this scaling factors to ad- 
just our numerical data onto one line for q = 1,2,3,4 
and L = 32,128. We plot / = \\og{it hs {p;L,q)) - 
\og(A(L, q))\ as a function of the new scaling variable 

Z = [B(L, q) (p - p c (L, q))f L - q) in FigE 

The scaling straight lines \z\ shown in Fig [3 corre- 
spond to the fitting function F(z) = exp(— |z|). We see 
that all points lie on the one curve, and this curve is very 
close to \z\ in the range 0.1 < z < 3. So, we might ex- 
pect that our fitting procedure is correct. However, on 
the graph, the tail points deviate from the line \z\, and 
this deviation explains the fact that moment ratios for 
TThsip', L.q) in the Table UTI differ from analytical values 
for g(x, C). Results of approximation for the crossing am- 
plitude A(L,q) are placed in Table ITTT1 For each value 
of q = 1,2,3 amplitudes A(L,q) are equal within our 
numerical accuracy of the approximation. Therefore, we 
can conclude that the scaling amplitude depends upon q 
and depends weakly upon L. 

Results of the approximation for the scaling index 
((L,q) are represented in the Table HVl We can see that 
this scaling index practically does not depend upon lat- 
tice size. We assume that for each value q we have the ad- 
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FIG. 2: The scaled crossing probability f(z) = | log(Q(z))| = 
| log(7Th s ) — \og(A(L, q))\ as a function of the new scaling vari- 
able z = [B{L,q){p - p c (L,q))] c{L ' q) for q = 1,2,3,4 and 
L = 32, 128. 



TABLE III: The scaling amplitude A(L, q) for 
16,32,64, 128 



iq=I;L=32 
1 q=l; L=128 
q=2; L= 32 
q=2; L=128 
q=3; L= 32 
q=3; L=128 
q=4; L= 32 
q=4; L=128 
Izl 



q 


1 


2 


3 


4 


A(L = 16) 


0.1806(2) 


0.166(1) 


0.151(2) 


0.137(2) 


A(L = 32) 


0.1794(3) 


0.167(1) 


0.151(2) 


0.132(2) 


A(L = 64) 


0.1806(3) 


0.167(1) 


0.148(1) 


0.128(2) 


A(L = 128) 


0.1794(3) 


0.166(1) 


0.149(2) 


0.122(2) 



ditional scaling index C{q)i which is nondependent upon 
the lattice size. Test values of £ in the first row of Table UTI 
are chosen to be close to ((q). It seems, that £(3) = 4/3, 
C(4) = 6/5. 

We know that in accordance with the scaling the- 
ory for each fixed value q the crossing probability must 

1 

be a universal function of variable x — (p — 

p c (L,q)). Therefore, we approximate the numerical data 
for B(L, q) (see FigEJ by function b{q)V^ and repre- 
sent results in Table Results of the approximation 



TABLE IV: The scaling index £{L, q) for L = 16, 32, 64, 128 



q 


1 


2 


3 


4 


C(L = 16) 


1.921(3) 


1.577(12) 


1.398(15) 


1.264(18) 


C(L = 32) 


1.900(4) 


1.557(12) 


1.364(13) 


1.238(14) 


C(L = 64) 


1.886(4) 


1.55(12) 


1.367(13) 


1.218(15) 


C(L = 128) 


1.887(5) 


1.545(14) 


1.337(13) 


1.198(14) 
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comparison. This fact once more confirms our approxi- 
mation procedure. The exception is the case q = 4, for 
which there is a difference between y(q — 4) and a . 

Many critical quantities in the Potts model q = 4 ex- 
hibit logarithmic corrections EE El EE El- These log- 
arithmic corrections explain the difference between ana- 
lytical value v (q = ±} = 1-5 and numerical approximation 
for the scaling index y = 1.39(6). 



IV. CONCLUSIONS 

To summarize results, we can conclude that the cross- 
ing probability TTh s (p; L, q) is a universal function Q(z) 



FIG. 3: Nominiversal metric factors B(L,q). Results of ap- 
proximation by the function b(q)L v ^ are shown by lines (see 
Table EJ. 



b(q)L v ^ are also shown in FigOUby lines. We can see 



TABLE V: The approximation of the nonuniversal metric fac- 
tors B(L, q) by the function b(q)L v(q) . The inverse correlation 
length index is added for the comparison. 



q 


1 


2 


3 


4 


b(q) 


1.737(6) 


1.65(3) 


1.51(3) 


1.38(3) 


y(i) 


0.756(8) 


1.011(3) 


1.218(16) 


1.39(6) 




0.75 


1 


1.2 


1.5 



that the thermal scaling index y(q) is very close to the 
analytical value of the inverse correlation length index 
jJ^y |23j , which is represented in the Table [V] for the 



iths{p;L,q) = A(q)Q f b{q)L~> 



&>(p-p c (L,q)) 



C(«) 



of the scaling variable z 



b(q)L^(p-p c (L,q)) 



(8) 



As we can conclude from Fig|2 the function Q(z) looks 
like the exponent Q(z) ~ exp(— \z\) on the interval 0.1 < 
z < 3, but deviates from it in the vicinity of and on the 
tails \z\ > 3. 

Let us pay attention to the new details of this work. 
We consider the universality of the crossing probability 
for different values of q by adding the scaling index £(<?). 
We work in the scale p, where p = 1 — cxp(— j3) is the 
probability of a bond to be closed instead of the usual 
scale t = (T — T c ) /T c to make the crossing probability 
symmetric. We find numerically the scaling index (,(q). 
The universal function looks like Q(z) ~ exp(— \z\) on 
the interval 0.1 < z < 3. 

The author would like to thank Robert M. Ziff for help- 
ful remarks, comments and discussion. The author is 
grateful to S. Nechaev for revision of the manuscript and 
to the Joint Supercomputer Center RAS (www.jscc.ru) 
for providing computational resources. 
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